Net quark number probability distribution near the chiral crossover transition 



CO 

o 

CN 

G 

c3 



Kenji Morita, 1, 2, B Bengt Friman, 3 Krzysztof Redlich, 4 ' 2 and Vladimir Skokov 5 

1 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
^Extreme Matter Institute EMMI, GSI, Planckstr. 1, D-6^291 Darmstadt, Germany 
'' GSI, Helmholzzentrum fur Schwenonenforschung, Planckstr. 1, D-64291 Darmstadt, Germany 
^Institute of Theoretical Physics, University of Wroclaw, PL-50204 Wroclaw, Poland 
5 Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 

(Dated: January 15, 2013) 

We investigate properties of the probability distribution of the net quark number near the chi- 
ral transition in the quark-meson model. The calculations are performed within the functional 
renormalization group method and under mean-field approximation. We show, that there is a sub- 
stantial influence of the remnant of the chiral phase transition on the properties of the probability 
distribution. In particular, at physical pion mass, the presence of mesonic fluctuations modifies 
the distribution from the mean-field one, as well as, from the Skellam distribution which is set as 
the reference for a non critical behavior. An intrinsic feature of the net quark number probability 
distribution is, that in the vicinity of the chiral crossover transition it is narrower than the corre- 
sponding Skellam function. We also discuss scaling properties of the probability distribution and 
the resulting cumulants. 
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I. INTRODUCTION 

The structure of the QCD phase diagram is one of the 
fundamental problems addressed in both, theoretical and 
experimental studies (HQ- At finite chemical potential, 
there have been predictions for the existence of the QCD 
critical point (CP), both within effective chiral models Q 
and some lattice QCD calculations Q . Despite extensive 
studies, however, the exact location of the CP and even 
its existence have not been established yet. 

So far, the lattice QCD simulation at physical quark 
masses has shown, that at vanishing baryon density the 
transition from hadrons to a quark matter is crossover, 
in which, no associated singularities exist in thermody- 
namic observables [f|. Therefore, possible experimen- 
tal signatures of the phase transition could be washed 
out. However, the chiral crossover transition in the QCD 
medium was conjectured [f| to appear as remnants of 
the second order phase transition belonging to the 0(4) 
universality class. Indeed, recent lattice QCD calcula- 
tions have confirmed the 0(4) scaling behavior of relevant 
physical observables at the physical quark masses @, [3 ■ 
This result has opened new opportunities to verify the 
QCD phase boundary experimentally by measuring fluc- 
tuations of conserved charges [9l-fl3j . 

The fluctuations of conserved charges, especially of the 
net baryon number and the electric charge, were pro- 
posed to search for the conjectured CP [9|, Il4| - tl8l ] . How- 
ever, in view of remnants of the 0(4) criticality and by 
virtue of the proximity of the chiral crossover to the 
freeze-out line in heavy ion collisions, fluctuations of the 
conserved charges were also proposed as diagnostic tools 
of the QCD phase boundary outside the CP [a El [II 01- 
l2ll . It has been shown, that at chiral crossover transition, 
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there should appear a non-trivial structure of the higher 
order cumulants of the net baryon number and electric 
charge due to 0(4) criticality [lo, HH- The theoretically 
expected modifications of different cumulants due to the 
chiral dynamics can be partly observed in recent data of 
STAR Collaboration [22|, [H| ■ The cumulants have also 
been studied in the first principle lattice QCD a ppr oach 
J2J-|26[, as well as, in effective chiral models [19W2U 

Fluctuations of conserved charges can be linked to the 
corresponding probability distribution. Thus, any criti- 
cal behavior in cumulants of conserved charges should be 
also reflected in their probability distributions. Recently, 
we have examined the effects of the phase transition on 
the net baryon number probability distribution based on 
the Laudau theory of the phase transition (34|. We have 
argued that there is a substantial influence of criticality 
on the probability distribution caused by the existence of 
the singularities of the thermodynamic potential in the 
complex chemical potential plane. In the present work, 
we extend our previous studies to more realistic model. 
We consider the two-flavor quark-meson model which is 
a low energy effective approach to the chiral properties 
of QCD. We include consistently the critical fluctuations, 
which were introduced ad hoc in Ref. [34[, by means of 
the functional renormalization group (FRG) method (35r - 
[37( | . Although we neglect coupling to the Polyakov loop, 
which is important to implement the statistical confine- 
ment, the 0(4) scaling behavior and thermodynamics 
near the critical point is well described by this model 
within the FRG method [H, H3 ■ Clearly, some proper- 
ties related with confinement in QCD, like e.g. relevant 
behavior of kurtosis near the chiral transition (l6H3l| can- 
not be addressed in this approach. Therefore, we restrict 
our studies entirely to effects of the chiral transition on 
the net quark number probability distribution. 

This paper is organized as follows: in the next section, 
we introduce the quark-meson model and its thcrmody- 
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namic properties. In Sec. IIII1 we present results on the 
probability distribution and cumulants of the net quark 
number. SectionQV]is devoted to the concluding remarks. 

II. THE THERMODYNAMIC POTENTIAL IN 
THE QUARK-MESON MODEL 

We employ the quark-meson model to explore the in- 
fluence of the chiral phase transition on the probability 
distribution of the net quark-number. 

The quark-meson model is an effective realization of 
the low energy QCD in which the meson field <f> = (a, tt), 
in the 0(4) representation, is coupled to quark fields 
through the covariant derivative. The Lagrangian den- 
sity is then given by 

C = q[i 7tl d» - g(a + i l5 f ■ 7r)]q + \{d^) 2 + \{d^f 
-U(a,Tf), (1) 



where U(a) denotes the mesonic potential, 

£/( ( 7,7r)-imV + ^ 4 -^. (2) 

For m 2 < and h = 0, the 0(4) symmetry of the po- 
tential is spontaneously broken to 0(3), resulting in a 
non- vanishing value of the vacuum scalar condensate (tr) . 
The last term, h = J^m 2 , breaks explicitly the chiral 
symmetry and yields the nonzero pion mass. 

To describe thermodynamics of the quark-meson 
model we calculate thermodynamic potential within the 
FRG approach as introduced in Ref. [38| . Following 
[35| , we consider a scale dependent effective action in the 
local potential approximation. We neglect the wave func- 
tion renormalization and the flow of the Yukawa coupling. 
Then, using the so-called optimized cutoff functions, the 
evolution equation for the scale dependent thermody- 
namic potential density with the reduced field varible 
p=(a 2 + tt 2 )/2 is found as H3, 



d k n k ( P ) 



12tt 2 



A {1 + 2n B (E„)} + i- {1 + 2n B (E a )} - ^ {l - 



E„ 



(3) 



where n B and uf are the Bose and the Fermi distribution functions, respectively, v q — 2N c Nf = 12 stands for the 
quark degeneracy The single particle energies of pion, sigma meson and quark/ antiquark read, 



fy> E a 



k2 + n' k + 2pti^ E± = yjk 2 +2g 2 p±p, 



(4) 



where f2J. and f2' fe ' denote the first and the second derivatives of Q k = ^ k + hy/2p kl with respect to p. The flow equation 
© is solved within the Taylor expansion method. Expanding the potential up to the third order in p around the 
potential minimum p kl 



a„(k) 



(P-Pk) n , 



(5) 



and using Eq. ((3j), one gets the flow equations for a n (k) 
and p k , which can be solved numerically starting at a cut- 
off scale A = 1.0 GeV [H|. We eliminate a\ by the scale 
independent relation h = ai(k)y / 2pk- Thus, there are 
four initial parameters in the flow equations to be fixed 
at the cutoff scale k = A. The ao(A) gives is a constant 
shift of f2. The finite momentum cutoff implies a non 
physical high temperature behavior of thermodynamic 
quantities due to the lack of hi gh momentum contribu- 
tion beyond the cutoff scale (3ll . |40| . Thus, following the 
prescription introduced in Ref. [40j, we include the high 
momentum contribution by considering the flow equation 
for the non-interacting massless quarks and gluons, 

d k Q.${T,n) = y^{2(^ c 2 - 1)[1 + 2n B {k)] 
- v„\l 



n F {k + )-n F {k~)]}. (6) 
Integrating the flow equation from k — oo to k = A, 



we obtain A (T, /z) which is used for ao(A). If one 
puts 0.3(A) = then the remaining two parameters, the 
running minimum p k =K arid 02(A) are fixed such that 
to reproduce the pion = 135 MeV and the sigma 
rn a = 640 MeV masses. The strength of the Yukawa 
coupling, g = 3.2, is fixed by the constituent quark mass, 
M q (T = 'p = 0) = g<J k = a {T = p = 0) = 300 MeV 
with a k=0 {T = p = 0) = U = 93 MeV. The full ther- 
modynamic potential density of the quark-meson model 
fi(T, /i), which includes thermal and quantum mesonic 
fluctuations, is obtained from tt(T,p) = limfc-^o^fcj 
where fl k =o is the solution of the flow equation (J3J). 

If one ignores the flow of the mesonic part in the evo- 
lution equation (|3]), the resulting effective potential cor- 
responds to that obtained in the mean-field approxima- 
tion with a finite cutoff. The fermionic vacuum fluctu- 
ation term, included in the mean-field (MF) potential, 
is necessary to reproduce the correct order of the phase 
transition at vanishing p. This vacuum term has to be 
renormalized to remove the artifact of the cutoff depen- 
dence [l||- Then, in the mean- field approximation, the 
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thermodynamic potential is obtained as 
n MF ({a);T, fi) = U((a),7f = 0) - ^L_M* In 



-u q T 



d 3 p 
(2^)3 



ln(l + e-^~^' T ) + ln(l + e^ E " + ^' T ) 

(7) 



where M q = g(a) and the expectation value (a) is deter- 
mined as the solution of the gap equation OCImf /d{d} = 
0. 



III. PROBABILITY DISTRIBUTION OF THE 
NET QUARK NUMBER 

The thermodynamical potentials for the quark-meson 
model, derived in the previous section, can be used to 
calculate the influence of the chiral phase transition, or 
its remnant, on different obscrvables. In the following we 
concentrate on properties of the probability distribution 
P(N) and the resulting cumulants c„(T, n) of the net 
quark number N . 



A. General properties of P{N) 

We consider a thermodynamic system described by the 
grand canonical ensemble at temperature T in a subvol- 
ume V. We introduce the chemical potential \i which is 
connected to the average net charge conservation. For 
the net quark number N = N q — Ng, the probability dis- 
tribution P(N) to find the net charge N in the volume 
V is given by 



P(N;T,^V) 



Z(T,V,N) /T 



(8) 



where Z(T, V, N) is the canonical and Z(T, fi, V) the 
grand-canonical partition function. Owing to the con- 
dition, X)jv=-oo -P(-^0 = 1' the grand canonical parti- 
tion function in Eq. ((5J) plays a role of the normaliza- 
tion. Thus, all essential information on P(N) proper- 
ties is contained in the canonical partition function. The 
Z(T, V, N) is connected with the grand canonical parti- 
tion function through the fugacity expansion, 



Z{T,V,\) = Y J ^ N Z(T,V,N), 



(9) 



N 



where A = e^ M . Consequently, the canonical partition 
function is obtained as the coefficient of this Laurent se- 
ries, 

z ^ v ^=^i dx ^W^- (io) 

Thus, to find the canonical partition function one needs 
to know the analytic structure of Z(T, V, A) in the com- 



plex A plane to choose an appropriate integration con- 
tour. 1 In the chiral effective models, the general structure 
of the singularities associated with the chiral phase tran- 
sition has been discussed in Refs. fill. |42|. In the broken 
phase, T < T c (fj,), there is no singularity on the unit cir- 
cle A = e ie with < 9 < 2ir. Consequently, the canonical 
partition function can be represented as [43l Eij , 



Z(T,V,N) 



1 

2^ 



271 



d6e- l6N Z(T, V, 



(11) 



where 8 = fij /T is the imaginary chemical potential. 

The above equation links the grand canonical parti- 
tion function at imaginary chemical potential to thermo- 
dynamics at a finite and fixed charge. The QCD par- 
tition function at imaginary chemical potential exhibits 
the Roberge- Weiss periodicity 0], Z(T, V 7 9 + 2tt/3) = 
Z(T,V,9), which implies the quantization of the net 
quark number by a multiple of three j46l - l48l |. Owing to 
the lack of confinement, the quark-meson model has the 
periodicity 2n in 9, thus here N is an i nteg er. 

In the present work, following Ref. [3J], we will use 
Eq. pip to calculate the canonical partition function and 
the corresponding probability distribution. The thermo- 
dynamic potential density 57 = —(T/V) logZ will be cal- 
culated in the quark-meson model within the FRG ap- 
proach and under the mean-field approximation. We also 
assume that the thermodynamic potential density is well 
approximated by its thermodynamic limit. Such an ap- 
proximation, also employed in Ref. [Hj], is justified for 
\N\ <C N*, where N* is a characteristic of the parti- 
tion function of a finite system in the Yang-Lee theory of 
phase transitions [iTjj . 

Because of an oscillating structure of the integrand, 
the numerical integration of Eq. f|lL[) has a limited ap- 
plicability at large \N\. We found, however, that the nu- 
merical integration works up to the value of \N\ at which 
P(N;fi = 0) ~ 10~ 12 , independent of the volume, tem- 
perature and other parameters. Such numerical accuracy 
is sufficient to study the influence of the chiral criticality 
on the properties of the net quark number probability 
distribution. 



B. The net quark number distribution near the 
chiral phase transition 

In the chiral limit and at vanishing and moderate val- 
ues of fj,, the quark-meson model exhibits a second order 
phase transition which belongs to the 0(4) universality 
class. For a physical pion mass, the chiral symmetry is 
explicitly broken. Nevertheless, there is a remnant of the 
0(4) criticality in different observables @. In particular, 



1 In a finite system, the grand partition function Z has Yang-Lee 
zeroes at complex A, which turn into cuts in the thermodynamic 
limit. 
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FIG. 1: Probability distributions of the net quark number just below the pseudocritical temperature T pc in the FRG approach 
(left) and in the mean-field (MF) approximation (right), compared with the corresponding Skellam distributions. The dots are 
the model results for different volumes while the lines stand for the Skellam distribution. 
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FIG. 2: Probability distribution for different temperatures in the FRG approach (left) and in the mean-field (MF) approximation 
(right). The dots are the model results while the lines stand for the corresponding Skellam distribution. 



such a remnant should be also present in the probability 
distribution of the net quark number distribution which 
is governed by the properties of the thermodynamic po- 
tential. Since the P(N) is calculated at finite volume, to 
see the influence of the long range critical correlations, 
the probability distribution should be calculated for a 
sufficiently large V and N and also near the pseudocrit- 
ical temperature T pc . 

To disentangle the influence of the chiral criticality 
on P(N), one needs to establish the reference. At low 
temperatures, T <C T pc , the thermodynamic potential is 
well approximated as a quasi-ideal quark gas with dy- 
namically generated mass which is T— dependent. Con- 
sequently, at any fixed T and V, the natural reference 
for P(N) is the Skellam distribution which reflects the 
probability distribution calculated from the non-singular 
part, an ideal gas approximation, of the thermodynamic 
potential . The Skellam distribution, in this case , 
is quantified by the mean number of quarks b = (N q ) and 



antiquarks b = (Nq) as, 

/,xJV/2 

P(N) = UJ I N {2Vbl>)e^ b+b \ (12) 

where In(x) is the modified Bcsscl function. 

In the quark-meson model, b and b are calculated for 
an ideal gas of constituent quarks with dynamically gen- 
erated mass 

v V f°° 

b= 2^J dfcfc2 M£fc;7», (13) 

where np is the Fermi distribution function and = 
■\Jk 2 + Mq 2 is the quasi-quarks energy. In the MF ap- 
proximation M q = g(cr), while in the FRG approach we 
apply the scale dependent quark mass M q ^ = guu which 
for k < A is obtained from the flow equation ([3]). For 
A < k < oo the quarks are assumed to be massless and 
non-interacting. 

Figure [1] shows the probability distributions of the net 
quark number obtained in the quark-meson model, both 
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FIG. 3: Probability distributions in FRG (red, circles), Skel- 
lam (blue, dashed), and free fermi gas (green, solid). 



in the FRG approach (left) and under the MF approxi- 
mation (right). The calculations were done at vanishing 
chemical potential in the vicinity of the pseudocritical 
temperature T pc . The values of T pc were determined by 
the peak position of the chiral susceptibility which in the 
FRG and MF cases correspond to 214 MeV and 190 MeV, 
respectively. P(N) is seen in Fig.[T]to be strongly chang- 
ing with the volume of the system. With increasing V 
the distribution is broader as expected, since the variance 
of the distribution is linearly dependent on the volume. 
There is also a strong change of the distribution when in- 
creasing distance from the pseudocritical point. As seen 
in Fig. [2] the distribution is narrower for smaller values 
of T/Tp C at a fixed volume. This is directly connected to 
a decreasing mean number of quarks in the system. 
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FIG. 4: Cumulants c„ (n = 2, 4, 6) of the net charge with the 
mean b = b — 6.0 calculated for the Skellam distribution from 
Eqs. (fT6)l - (fT9)l for different values of iV max . 

In comparison, both the MF and the FRG distribu- 
tions are narrower than the Skellam counterparts. The 



deviation becomes larger when the volume increases and 
the mesonic fluctuations are included, as seen in Fig. [1] 
The effect of the mesonic fluctuations is more clearly vis- 
ible when temperature is varied. In Fig. [21 we plot P{N) 
for three different temperatures below T pc . Except for 
the highest temperature, the P(N) in the MF calcula- 
tions coincides with the Skellam distribution, while it 
deviates in the FRG case. As temperature approaches to 
Tp C , the difference between the Skellam and the model 
results becomes larger, indicating an increasing influence 
of remnants of the 0(4) criticality due to the chiral dy- 
namics. 

While the shrinking of the P(N) in the FRG cal- 
culations from the Skellam distribution could be ex- 
pected from the previous studies of the Laudau theory 
in Ref . [34| , one may ask why such behavior is also seen 
in the present MF approach. As discussed in j34[, the 
shrinking of P(N) originates from the branch point in a 
complex chemical potential plane and its associated cut, 
near which, the two saddle points are located. Although 
the <7 4 model used in Rcf. [HJ does not have a cut in 
thermodynamic potential but in the order parameter, in 
general, however, such a cut exists also in the thermo- 
dynamic potential under the MF approximation of the 
quark-meson model, leading to the observed shrinking of 
P(N) below the Skellam distribution. 

The Skellam distribution (fT2"|) . corresponds to a non- 
interacting gas of classical particles. In general, effects 
of quantum statistics should be taken into account when 
the mass of quark is small relative to the temperature. In 
Eq. ([T3]), we use the Fermi distribution function, thus the 
effect of quantum statistics is partly taken into account 
through the number of quarks and antiquarks. Never- 
theless, the quantum statistics can also change the func- 
tional form of the probability distribution, as shown for 
the pionic electric charge fluctuations in Ref. [ll| . In the 
quark-meson model within the FRG calculation one gets, 
M q ~ 126 MeV at T/T pc = 0.98, the mass, which is too 
small to verify the Boltzmann approximation near T pc . 
Therefore, to disentangle the criticality and the quantum 
statistics effect, we also calculate the probability distribu- 
tion P(N) for the free quark and anti-quark gas through 
Eq. (jlip from the thermodynamic potential density 



+ \n{l + e- E "l T - i6 )] 



-E k /T+ie 



(14) 

with the same E}. as used in Eq. (fT5|) . Note, that by 
taking the Boltzmann approximation in Eq. ()14[) one re- 
covers the Skellam distribution (|12[) in which the baryon 
number is given by the Boltzmann distribution. 

Figure [3] shows, that the probability distributions 
P{N) for the free Fermi gas exhibits narrower width than 
the corresponding Skellam distribution, reflecting the in- 
fluence of the quantum statistics. However, the P{N) 
obtained in the quark-meson model within the FRG ap- 
proach is still narrower than for the corresponding Fermi 
gas, indicating the effect of the 0(4) criticality. 
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C. Cumulants of the net quark number 



Cumulants, c n (T,/j,) of the net-quark number, arc 
known to be excellent probe of the chiral phase transi- 
tion and its remnant d, [HI [H, [H-H3 • In the statistical 
physics, c„ are quantified by the corresponding suscepti- 
bilities, 



d n {p(T,Li)/T 4 
dQi/T)» 



(15) 



Thus, once the thermodynamic pressure p(T, /i, V) = 
(T/V) In Z is calculated in the grand-canonical ensemble, 
the cumulants c„ (T, fx) can be obtained by taking deriva- 
tives of the thermodynamic pressure with respect to the 
chemical potential. Clearly, since the higher order cumu- 
lants are increasingly sensitive to the chiral phase tran- 
sition, their calculations involve special numerical efforts 
[1, H3, E3l- However, once the probability distribution of 
the net charge multiplicity, Eq. (|8j , is known, the cumu- 
lants can be calculated directly from P(N) in terms of 
powers of the central moments ((SN) k ) — ((N — (N)) k ). 
Indeed, at \x = 0, the first three non- vanishing cumulants 



c„ (n = 2, 4, 6) are expressed by the central moments as 

<(<w) 2 ) 



C2 



VT 3 ' 
((SN) 4 ) - 3((<WV) 2 ) 



(16) 
(17) 



VT 3 

c 6 =[((6N) 6 )-15((5N)*)((5Nf) 

-10({SN) 3 ) + 30({5N) 2 )] /{VT 3 ); (18) 

where (N k ) are directly obtained from the probability 
distribution 



(N k ) 



N k P(N), 



(19) 



N=-N„ 



where in general, iV max = oo. In practice, however, the 
N max is always finite. Consequently, the cumulants cal- 
culated from Eqs. (|Tfj|) - (fT5|) are only approximations of 
exact results obtained by taking derivatives of the ther- 
modynamic pressure as in Eq. (|15|) . To illustrate this 
property we calculate different cumulants of Skellem dis- 
tribution for different iV max and compare those results 
with exact calculations. Figure2]shows the ratio of differ- 
ent moments as a function of N max for a fixed mean num- 
ber of quarks. It is clear from this figure, that to repro- 
duce the exact results for different cumulants, one needs 
to calculated P(N) for sufficiently large N max , which 
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is increasing with the order of the cumulant. Clearly, 
the value of iV max which is required to get the best de- 
scription of different cumulants also varies with the mean 
charge number. In the actual model calculations, the 
mean quark number depends, in general, on the volume, 
the temperature and the chemical potential. 

Figure [5] shows the temperature dependence of dif- 
ferent cumulants calculated in the quark-meson model 
within the FRG approach from Eq. (|15|) and their approx- 
imate values obtained from Eqs. ([T6f -([T9 f for different 
N max . The cumulants obtained from the correspond- 
ing Skellam distribution are also shown in this figure for 
comparison. 

The convergence properties of different cumulants to 
their exact values with iV max , in the quark-meson model, 
are similar to that already seen in the Skellam distri- 
bution. The A max needed to get a good approximation 
of the exact results are increasing with the order of cu- 
mulants, since they are more sensitive to the tail of the 
distribution. For the parameters used in Fig. [5j all cu- 
mulants, up to the sixth order, are well reconstructed 
from the probability distribution by choosing -/V max ~ 90. 
This also confirms consistency of the calculations of the 
canonical partition function. 

At fj, = 0, the second and the fourth order cumulants 
are not influenced by the critical chiral dynamics, as they 
are finite even in the chiral limit. Their properties are 
entirely determined by the regular part of the partition 
function. The temperature dependence of C2 and C4, seen 
in Fig. [5J is due to the temperature dependence of the 
dynamical quark mass. 

Contrary to c-i and C4, the temperature dependence of 
C6 and its negative values near T pc , seen in Fig. [5j has 
a universal structure. The characteristic shape of cq is 
typical for the 0(4) universality class, originated from the 
universal scaling function [i| . The cq obtained from the 
noncritical Skellam distribution, exhibits a very different 
behavior. 

It is interesting to note, that already at moderate val- 
ues of -/V max , the 0(4) shape structure of cq can be re- 
produced from the probability distribution. This is an 
interesting result for heavy ion experiments where one 
expects to observe the chiral criticality by reconstructing 
the higher order moments from a measured probability 
distribution of the net charge [22| . 

Deviations of c^, obtained from the probability dis- 
tribution (fT6]) - ([T9")) , from their exact values c„, obtained 
from Eq. ([5]), also depend on the volume of the system. 
However, there is an interesting approximate scaling re- 
lation (APS) which describes the relative deviation of 
from c^f for different volumes. Figure [5] shows the relative 
difference, R = (c^ — c^)/c^ for n = 2, 4, 6, as a function 
of A max normalized by y/V. There is a clear APS of c-2 
and those of higher order cumulants for sufficiently large 
A max / vF. Since for fixed T the mean number of quarks 
N q scales linearly with the volume, one can replace y/V 
by y/N q to preserve the observed approximate scaling 
properties. 
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FIG. 7: Rescaling of the net quark number probability dis- 
tribution at fixed T with volume in the quark-meson model 
within the FRG approach and for the corresponding Skellam 
distribution. 



The moments of the net quark number are linked to the 
corresponding probability distribution through Eq. (|I9I) . 
Thus, the APS of cumulants, seen in Fig. |Hl should be 
also reflected in the properties of the net quark number 
probability distribution. Figure [7] shows, that for fixed 
T the rescaled probability distribution y/V P(N) indeed 
scales approximately with N/y/V. The APS is valid 
for P(N) obtained in the quark-meson model within the 
FRG approach, as well as, for the corresponding Skellam 
distribution. Therefore, the strong variation of P{N) 
with V, seen in Fig. [TJ can be essentially reduced if plot- 
ted as functions of the rescaled variable, N/y/V. The 
results in Fig. [7] also confirm, that shrinking of the prob- 
ability distribution of conserved charge, relative to the 
Skellam distribution, due to the 0(4) chiral dynamics is 
almost independent of the system volume. 

Figure [5] shows probability distributions rescaled 
with the standard deviation of the distribution, a = 
y/{(SN) 2 ) = y/c 2 VT i . Probability distributions ob- 
tained from the model calculation and the correspond- 
ing Skellam functions are hardly distinguishable on the 
logarithmic plot. They are also very different from the 
Gaussian distribution. The above results imply, that the 
global shape of the probability distribution with 0(4) 
criticality is not significantly different from the Skellam 
distribution, in spite of a very different behavior of the 
higher order cumulants. However, when approaching to 
the critical point, there is a clear deviation from the Skel- 
lam distribution even after rescaling. 

Figure [5] demonstrates the violation of the "scaling be- 
havior" close to the second order phase transition, in 
which the pion mass was reduced according to h = 
10 _8 / w m^.. There is clear anomalous behavior in the 
tail of the distribution, which can be interpreted as the 
realization of the oscillation due to the saddle points dis- 
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FIG. 8: Rescaling of the probability distribution with the 
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the FRG distribution in the quark-meson model and for the 
corresponding Skellam distribution. The standard normal dis- 
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line. 
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of the critical chiral dynamics governed by the 0(4) uni- 
versality class in terms of the quark-meson model. Our 
calculations of P(N) were done within the Functional 
Renormalization Group (FRG) approach which preserves 
the 0(4) scaling of relevant physical observables. 

The main objectives of this paper was to verify the 
influence of the remnant of the chiral phase transition, 
expected at the physical pion mass, on the net quark 
probability distribution. We have shown that P(N) is a 
sensitive observable to criticality by comparing the FRG 
results to the mean-field (MF) calculations. The thermal 
and quantum mesonic fluctuations, which are present in 
the FRG approach, imply modifications to the distribu- 
tion obtained under the MF dynamics. 

We have shown, that near the chiral crossover transi- 
tion, the probability distribution of the net quark num- 
ber is narrower than the Skellam distribution, which was 
used as the reference for the non-critical distribution at 
the corresponding temperature. Such narrowing of the 
probability distribution is dominated by the 0(4) criti- 
cality. In the quark meson model, however, it can also ap- 
pear due to the quantum statistic effect as a consequence 
of a decreasing quark mass near the chiral transition. In 
the QCD medium, the masses of baryons are not directly 
effected by the critical chiral dynamics. Therefore, the 
narrowing of the net baryon number or the net proton 
number probability distribution measured in heavy ion 
collisions can be directly linked to the 0(4) criticality. 

We have found an approximate scaling properties of 
the probability distribution and the net charge cumu- 
lants for different variables and thermal parameters. We 
have shown, that there is approximate scaling of P{N) 
with the volume of the system, implying that the ob- 
served shrinking of P(N) near the chiral transition is 
very weekly volume dependent. We expect, that our re- 
sults are of importance in heavy ion collisions, where by 
measuring the net baryon number probability distribu- 
tion and related moments, one expects to probe the QCD 
phase boundary experimentally. The phenomenological 
implications of our results will be presented elsewhere 



FIG. 9: Probability distribution of the quark-meson model 
in the chiral limit rescaled with the standard deviation for 
different volumes. 
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